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I. INTRODUCTION 

In recent years the search for nonintegrabhhty cri- 
teria for Hamiltonian systems in the complex domain has 
acquired more relevance [lj]-[5j. Such techniques are po- 
tentially of particular importance in cosmology because 
of controversies over both integrability, and the existence 
of chaos in cosmological models [6]-[14]. Part of the prob- 
lem is that certain methods used to traditionally mea- 
sure chaos in non-relativistic systems, such as the Lya- 
punov exponents, are no longer valid in General Relativ- 
ity where there is no absolute time coordinate. 

In this work we use a recently developed theorem by 
Morales and Ramis [l[ which establishes a relation be- 
tween two differents concepts of integrability: the com- 
plete integrability of complex analytical Hamiltonian sys- 
tems (given by Liouville's theorem) and the integrabil- 
ity of homogeneous linear ordinary differential equations 
(LODEs) in terms of Liouvillian functions in the com- 
plex plane. A Liouvillian function is a function which 
can be written as a combination of elementary func- 
tions, algebraic functions (solutions of polynomial equa- 
tions), their indefinite integrals or exponentials of these 
integrals. Since we are working in the complex domain 
this definition includes elementary functions such as log- 
arithms, trigonometric functions and their inverses. 

The model we study in this paper is a Friedman- 
Robertson- Walker (FRW) cosmological model with a 
conformally coupled scalar field, (f>, of mass m. We use 
the conformal form of the metric 



1 - kr 



■ dr^ - dfl^ 



(1) 



with Tj the conformal time, a{ri) the scale factor and k = 
0, ±1 the curvature. 

The dynamics of this model has been discussed and 
studied before using numerical methods but, as 

far as we are aware, no completely rigorous conclusion 
has been reached about its integrability. In the in- 
tegrability of a generalisation of the model studied here 
is considered using Painleve analysis via the ARS algo- 
rithm Though there is a strong connection 
between integrability and the Painleve property, and the 
latter has been remarkably successful in indicating possi- 
bly integral cases, it is worth noting that the lack of the 
Painleve property is not a rigorous obstruction to inte- 
grability Additionally, the ARS algorithm is not a 
foolproof method for determining whether a system pos- 
sesses the Painleve property, and its application can lead 
to false conclusions pH [Tsj . particularly when ap- 
plied to determining the non- integrability of a dynamical 
system. Finally, certain expressions in are undefined 
for our model, requiring a separate analysis. 

The Morales-Ramis Theorem (MRT) which we use in 
our study rigorously provides necessary conditions for 
the integrability of a Hamiltonian system and so suffi- 
cient conditions for non-integrability. The theorem can 
be used to reduce the question of integrability to one of 
the existence of Liouvillian solutions of a homogeneous 
second-order linear ODE. This problem can in turn be 
solved using Kovacic's algorithm. Though complex to 
write down, the algorithm is, as we shall see, straightfor- 
ward to apply to the problem considered here. To use 
the MRT for our cosmological model we first require a 
Hamiltonian which generates the field equations. In this 
case it is known |6l-|9| that a suitable Hamiltonian is 
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[{pI + kcl,^) - {pI + ka^) + m^a^c^^] = 0, (2) 

where Pa and p^ are the momenta conjugate to a and 
respectively. 
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In the next section we give the main results of the 
Morales-Ramis theorem. Since the various versions of 
Kovacic's algorithm in the literature [H, [3, [13, HH have 
slight differences in presentation and conventions, we in- 
clude the version of the algorithm as used by us. We 
then show how the algorithm quickly determines that 
the Hamiltonian is nonintegrable for fc 7^ 0. Finally for 
the case k = the analysis based on the invariant planes 
a ^ Pa — and 4> — p,f, — is inconclusive. However, 
because the potential is homogeneous in this case, there 
exist particular nonsingular solutions which do not lie in 
these planes which can be used as a basis for the analy- 
sis. Fortunately the case of homogeneous potentials has 
been exhaustively studied by Yoshida Q and Morales- 
Ramis [T] and so we can simply apply those results. 



II. THE MORALES-RAMIS THEOREM 

The Morales-Ramis theorem is a nonitegrability cri- 
terion: it gives a necessary condition for a Hamiltonian 
system to be integrable and therefore a sufficient condi- 
tion for nonintegrability. The theorem is based on the 
analysis of the variational equations (in particular the 
normal variational equation, or NVE) for the perturba- 
tions of a non-equilibrum particular solution. The basic 
idea is that if the flow of the Hamiltonian system has 
a regular behaviour (is integrable), then the linearized 
flow along a particular integral curve given by the NVE 
must also be regular (integrable). Conversely if the lin- 
earized flow is nonintegrable the system as a whole will 
be nonintegrable. 

A Hamiltonian system, Xh, of dimension n is called 
integrable if there exist n independent constants of the 
motion in involution. By considering the differential Ga- 
lois group of the NVE, the theorem of Morales-Ramis 
links this concept of integrability to an apparently differ- 
ent concept of integrability - the existence of Liouvillian 
solutions of the NVE of Xh- The theorem may be stated 



Theorem 1 // there are n first integrals of Xh that are 
independent and in involution, then the identity compo- 
nent of the Galois group of the NVE is abelian. 

It is known that 19] for an ODE to admit a Liouvillian 
solution, the identity component of its Galois group must 
be soluble. Hence, if the solutions are not Liouvillian, 
the identity component of the Galois group is not soluble 
and, therefore, non-Abelian. 

Our strategy will therefore be: 
1: Select a particular solution (in our case an invariant 
plane) . 

2: Write the variational equations and the NVE. 
3: Check if the solutions of the NVE are Liouvillian func- 
tions. 

To decide the third step, we use Kovacic's algo- 
rithm }2u\ which we now turn to describe. 



III. KOVACIC'S ALGORITHM 

Kovacic's algorithm provides a procedure for com- 
puting the Liouvillian solutions of a homogeneous linear 
second order differential equation. If the algorithm termi- 
nates negatively, we can conclude that no such solutions 
exist. 

Let C(a;) be the field of rational complex functions (ra- 
tios of polynomials in x with complex coefficients). It is 
well-known that by using the change of dependent vari- 
able 



y = C exp 



5 da; 



(3) 



the second order homogeneous LODE 

y" + b{x) y' + c{x) y = 

can be transformed to the so-called reduced invariant 
form 



where 



e" - 5^ = 0, 



six) = h'ix) + h{x)^ 



c{x). 



(4) 



(5) 



Note that, if b{x) and c{x) S C(a;) then g{x) £ C{x). 

Moreover, using a further change of variables v = ^' /£,, 
equation ([4]) is transformed into the Riccati equation 

V + v'^ — g. (6) 

Now equation ^ is integrable, if and only if equa- 
tion ([6]) has an algebraic solution, that is v solves a 
polynomial equation f{v) = 0, where the degree of / 
(the minimal polynomial) in v belongs to the set L = 
{1,2,4,6,12}. 

Kovacic's algorithm can be divided into three main 
steps: the first step is the determination of the subset 
of L relevant for the LODE under consideration; the two 
other steps are devoted respectively to determining the 
existence of the minimal polynomial, and its construc- 
tion. If the algorithm does not terminate successfully (ie, 
equation ([6]) has no algebraic solution) then equation (|4|) 
has no solution in terms of Liouvillian functions. 

In the version used of the algorithm we essentially fol- 
low [EQUHIll. Let 
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t{x)' 



(7) 



with s{x),t(x) relatively prime polynomials, and t{x) 
monic. Define the function h on the set Lmax — 
{1,2,4,6,12} by hil) = 1, hi2) = 4, /^(4) = /i(6) = 
h{12) = 12. 

Step 1 (determination of possible orders of the mini- 
mal polynomial) 
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If t{x) = 1 then set m = 0, else factorize t{x) into 
monic relatively prime polynomials 

where U have no multiple roots and 7^ 1- 
Then 

1.1 Let r' be the set of roots of t{x) (i.e., the singular 
points in the finite complex plane) and let F = F' U 00 be 
the set of singular points. 

Then the order of a singular point c e F' is, as usual, 
o(c) = i if c is a root of multiplicity i of t^. The order at 
infinity is defined by o{oo) = max(0, 4 + deg(s) — deg(t)). 
We call m+ = max(m, 0(00)). 

For ^ i ^ m+, denote by F^ = {c G F | o(c) = i} the 
subset of all elements of order i. 

1.2 If m+ > 2 then we write 72 — card(F2), else 72 = 0. 
Then we compute 



2.5.1 If c G F', 



gc 



7 = 72 + card 



u 



Ffc 



/ 

1 

0(1), 

Oix~^) 



3 < fc < m+ 

^ k odd 

1.3 For the singular points of order one or two, c G 
F2 U Fi, we compute the principal parts of g: 

gc = ac{x - c)"^ + I3c{x ~ c) 
if c G F', and 

goo — (^coX 4~ I^oqX 

for the point at infinity. 

1.4 We define the subset L' (of all possible values for 
the degree of minimal polynomial) as {1} C £' if 7 = 72, 
{2} C L' if 7 > 2 and {4, 6, 12} C L' if m+ < 2. 

1.5 We have the three following mutually exclusive 
cases: 

1.5.1 If m+ > 2, then £ ^ L'. 

1.5.2 Define = y^l + Aa^. If m+ < 2 and Vc G 
Fi U F2, Ac G Q, then L = L'. 

1.5.3 If cases (1.5.1) and (1.5.2) do not hold, then L = 
L'- {4,6,12}. 

1.6 If L — 0, then equation ([4]) is non-integrable with 
Galois group SL(2, C), else one writes n for the minimum 
value in L. 

For the second and third steps of the algorithm we 
consider a fixed value of n. 
Step 2 

2.1 If 00 has order we write the set 

h{n)' 



n n n n 

2.2 If c has order 1, then Ec — {h{n)}. 

2.3 If n = 1, for each c of order 2 we define 

i?c = |i(l + A,),i(l-A,) 

2.4 If rt > 2, for each c of order 2, we define 



= z n ^ ^(1 - A,) + ^kAc : fc - 0, 1, . 



2.5 If = 1, for each singular point of even order 2i>, 
with v > 1, we compute the numbers ac and (3c defined 
(up to a sign) by the following conditions: 



[ (a; - c) 
and we write 



i=2 



7-^1 + — ^^ + oix-cy 

{x-cy\ (a;-c)"+' 



^/g^ := ac{x - c) " + ^ Mi,c 
2.5.2 If c = oo, 



X — c 



i=2 



9oc — i^Q^oo^ 

and we write 



u-2 




p^x--^ + O {x''-'') 



'9o 



l-^i, ooX 



i=0 



Then for each c as above, we compute 
and the sign function on Ec is defined by 

being +1 if = 0. 

2.6 If n = 2, for each c of order ly, with > 3, we write 
Ec = M. 

Step 3 

3.1 For n fixed, we try to obtain elements e — (ec)cer 
in the Cartesian product Hcsr ^c, such that: 

(i) d(e) := n — ^ Cc is a non-negative integer. 



h{n) 



cer 



(ii) If n = 2 or n = 6 then e has an even number of 
elements which are odd integers 

(iii) when n — A, then e has at least two elements not 
divisible by 3, and the sum of all elements not divisible 
by 3 is divisible by 3. 

If no such set e is obtained, we select the next value 
in L and repeat Step 2, else n is the maximum value in 
L and the Galois group is SL{2, C) (and equation ^ is 
non-integrable) . 

3.2 For each family e as above, we try to obtain a 
rational function Q and a polynomial P, such that 

(i) Q = TTTT X! 7~r7 + '^"1 X! sign(ec)v^. 



h(n) — ' X — c 



where (5„i is the Kronecker delta. 

(ii) P is a polynomial of degree d(e) and its coefficients 
are found as a solution of the (in general, overdetermined) 
system of equations 

P-i = 0, 

- -(P,)' - QP, -{n-^){i^- l)gP^+i, n > i > 0, 
Pn = ~P- 

If a pair (P, Q) as above is found, then equation ^ is 
integrable and the Riccati equation ([6]) has an algebraic 
solution V given by any root v of the equation 

n 

^-^ in — ly. 
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If no pair as above is found we take the next value in 
L and we go to Step 2. If n is the greatest value in L 
then the Galois group of (g]) is 5-^(2, C) and the ODE is 
non-integrable. 



IV. APPLICATION AND RESULT 

A. The case fc / 

We apply the theorem of Morales-Ramis to ([2]). We 
choose as our set of non-equilibrium particular solutions 
the invariant plane pa = a = 0. The NVEs relative to 
this plane are 

^ = H + (8) 

Changing the independent variable to (j) and renaming 
5a — y, we obtain the equation 

d^y Idy (m^ 1 \ 

This equation is a second-order, linear and homoge- 
neous ODE with coefficients which are rational functions 
of 0, and we can therefore apply Kovacic's algorithm to 
determine any Liouvillian solutions. 

Using ([3]) and ^ we transform ^ into the reduced 
invariant form ([4]) 



Summarizing, 

i?o = ||-^| and £;oo = {l}. 
Step 3 

In this step we need to calculate HcGr -^c, using the 
sets determined in the previous step. We obtain the set 
of sets given by 




From 3.1(i) we calculate the values of d[<s) as = — - 

and — 2 '^^spectively Since neither of theses values sat- 
isfies 3.1(i) and there are no other values of n in L, the 
Galois group of (fTO| is SL{2,C), equation pO| is nonin- 
tegrable in terms of Liouvillian functions, and therefore 
the system represented by the Hamiltonian ([2]) is also 
nonintegrable when k ^ 0. This completes the proof. 



B. The case fc = 



Lemma 1 Equation (0) has no Liouvillian solutions 
when k ^ and ni ^ 0. 

Proof: by application of Kovacic's algorithm to equa- 
tion uni). 

Step 1 

has one finite pole, at </) = 0, of order 2 and 
the pole at infinity, of order 4 (since, by assumption, we 
are treating the massive case, to 7^ 0). This implies that 
= 4 and 7 = 72 = 1 • Since the pole at = belongs 
to we calculate the Laurent series (when fc 7^ 0) as 




Hence ao = - and /3o — 0. Thus we have L — {1}. 
Step 2 

Because L = {1} the unique value for n n — 1. 
Through the items 2.3 and 2.5 we calculate the sets i?c. 

From 2.3 we have that En — l"^, — 1. In item 2.5.2 we 

\2 2/ 

need to expand g around (f) = 00. Doing this we obtain 



Lemma 2 When fc = the only first integral of the 
Hamiltonian system (0j is the Hamiltonian, and the sys- 
tem is therefore nonintegrable. 

Proof: In order to prove this lemma we observe that 
when fc = the Hamiltonian Q has a homogeneous po- 
tential. Hamiltonians with homogeneous potentials have 
been exhaustively studied using the MRT and particular 
results obtained which we now outline. 

Let 



H=-Y,Pl + V{qi,...,qn) (11) 

1=1 

where A is a constant and V is a homogeneous potential, 
i.e. V{AQ) = A^ViQ) with g being the degree of the 
potential. To put our Hamiltonian in the form (|lip we 
perform the canonical transformation, x = ia and Px — 
—iPa, after which 

H=]^[pl+pl-m^x^''] (12) 

and = 4. The MRT for homogeneous potentials is given 
by 
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Theorem 2 Let V{qi, . . . , qn) be a homogeneous poten- 
tial function of integer degree g, c a solution of the equa- 
tion c = V'{c), and Xi (the Yoshida coefficients) the 
eigenvalues of the matriz V"{c). One of these eigenvalues 
is trivial, in that it corresponds to the tangential varia- 
tional equation, and has value g — 1- 

If a Hamiltonian system of the form ill]) is completely 
integrable (with holomorphic or meromorphic first inte- 
grals) then each pair {g, Xi) belongs to one of the following 
list ( where we do not consider the trivial case g — 0) 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 

(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



g,P + p{p- 1)5/2) 

2, arbitrary complex number) 

—2, arbitrary complex number) 

+ 10p)2) 

4p)2) 
+ 6p)2) 
+ 6p)2) 
+ 6p)2^ 



'10 

-5^-^(4 

'8 8 V 



24 

-3, 



8^3 ' 

M2- 



24 

24 ^ 

— — + 

24 ^ 

— — + 

24 ^ 

— — -I- 

24 ^ 

_ii 

40 T 
2V g 



24 

J_C3 

24 V 5 
2^(2 

-(- 



-6p)2) 
6p)2) 
6p)2^ 



(13) 



24 



J_(10 
4,0 V 3 

3b(4 + 10p)2 
+ P(P 



10p)2) 



1)5)) 



where p is an arbitrary integer. 



For the system represented by ()12|) the only non-trivial 



Yoshida coefficient is A = —1. For 5 = 4 the only pos- 
sibilities for satisfying Theorem 3 are (1), (15) and (18). 
For all these cases there are no integer values of p which 
solve A = — 1, and we conclude that the system repre- 
sented by (|12p is nonintegrable. This completes the proof 
of lemma 2. 

We can now enunciate the following theorem 

Theorem 3 The Friedmann Robertson Walker model 
with a conformally coupled massive scalar field repre- 
sented by the Hamiltonian (0) is not completely inte- 
grable. 

Proof: By lemma 1 there are no Liouvillian solutions of 
the NVE for the plane a = Pa = 0. This implies that the 
identity component of its Galois group is not soluble and 
therefore non-Abelian. Using theorem 1 we have that 
the only first integral of the Hamiltonian system is the 
Hamiltonian itself, and that the system is not completely 
integrable in this case. By lemma 2 the Hamiltonian 
system is also not completely integrable when fc = 0. 
Therefore the Hamiltonian system represented by ([2]) is 
nonintegrable for all values of fc. 

V. CONCLUSION 



From our analysis we have shown rigorously using an- 
alytic methods that FRW universes with a conformally 
coupled massive scalar field are nonintegrable. This is 
compatible with results from numerical analysis based on 
Poincare sections [61] which indicate that the behaviour 
of the system is mathematically chaotic. 
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